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SUMMARY 


This paper is concerned with the determination of the numerical values of some of 
the eigenvalues and the total number of eigenvalues of the differential equation x 2 p” + 
[2dxe -x - k (k+1) - u 2 d 2 x 2 ] P = O. The problem arises in treating the effective ioni- 
zation potential of an atom in the interior of a plasma, but is essentially a mathematical 
problem. 

The problem is first treated from a theoretical standpoint. Then numerical tech- 
niques are developed that lead to approximations of the eigenvalues. An approximate 
inequality is obtained for the number of eigenvalues and an approximate equation is ob- 
tained for the eigenvalues. Tabular data are shown listing a total of 55 of these eigen- 
values for nine different cases. 


I. INTRODUCTION 

The problem treated in this note arises in the context of theoretical physics. It 
is concerned with the partition function of an atom in the interior of a plasma as dis- 
cussed by G. Ecker and W. Weizel. (See the Reference at the end of this note. ) The 
following is quoted from the cited paper: "As regards the partition function there is a 
well-known divergence difficulty whose elimination is of interest both numerically and 
as a matter of principle. Existing attempts to limit the number of terms either start 
from impermissible assumptions or content themselves with a crude cutting-off of the 
series of terms. Measurements of the effective ionization potential by various proce- 
dures have given results that are in marked contradiction with each other. In the way 
of a theoretical determination of the lowering of the ionization potential, only rough es- 
timates are available. " The difficulty mentioned in this quotation is primarily that of 
determining the eigenvalues of the differential equation 
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where m, h, E, Z, D, L are constants having special physical significance as explained 
in the cited paper. 

This note will be restricted to a mathematical treatment of the problem of deter- 
mining the eigenvalues of a differential equation equivalent to the one above. The methods 
used are applications of classical procedures commonly used in mathematics. Unless 
explicitly defined, no symbols will be used other than the usual symbols of mathematics. 

Acknowledgement is made to Mr. Neal Shepard and to Mr. A. Anderson of the 
General Electric Computation Unit at Marshall Space Flight Center in Huntsville, Ala- 
bama for their work in programming the evaluation of large determinants and a routine 
based on a recursion formula, respectively. 


The differential equation mentioned in the introduction can be transformed by sub- 
stitution to the following form: 


where d, k, u are constants and the equation is subject to the boundary conditions P(x) 

— O as x —O and P(x) — O as x — °°. The parameter d is to assume positive values 
which may be taken as integral powers of ten; the parameter k is to assume non-negative 
integral values; and the parameter u 2 gives the eigenvalues that will cause the solution 
P(x) to satisfy the boundary conditions. It may be noted that a sign factor has been in- 
troduced so that the eigenvalues will be positive in our treatment, whereas they are 
negative in the treatment by Ecker and Weizel, which treatment concerns the case k — O 
and eliminates the term involving k in equation i. . 

This equation is a second order differential equation with a regular singular point 
at the origin and no other singularities in the finite complex plane. There is an irreg- 
ular singular point at infinity and the object is to find the values of u 2 (or equivalently 
the values of u) such that the boundary condition P(x) —O as x — °° along the real axis 
is satisfied. The solution P(x) = O is ruled out as trivial in the following. 

The roots of the indicial equation are readily found to be -k and k + 1. Since these 
are integers, only the largest, namely k + 1, leads to an analytic solution. Thus the 
smallest exponent of x in the power series for P(x) must be k + 1. From this it follows 
that P(O) = O as required by one of the two boundary conditions (since k = O, 1, . . . ) . 


H. DISCUSSION OF THE DIFFERENTIAL EQUATION 


x 2 P" + [ 2dxe X - k (k+1) - u^x 2 ] P = O, 


(1) 
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From equation 1 and the condition P — O as x — °°, we see that P" - u 2 d 2 P = O 
will roughly determine P for extremely large values of x. Since the general solution of 
this equation is ae llC ^ x + be _u< ^ x , the condition that is to be satisfied requires the choice 
of e -uc ^ x as a rougn approximation to P(x) near infinity. This suggests that the substi- 
tution P(x) = e -u dxy be m ade in equation 1. This leads in a straightforward manner to 
the equation 

x 2 y" - 2 udx 2 y’ + [2dxe -x -k (k+1) ] y = O. (2) 


Taking into account the substitution and our previous knowledge about equation 1 , we sub- 
stitute the series y = T^aj x^ + ^into equation 2, which gives the following recursion for- 
mula for the determination of the coefficients in the power series for y. 


i (i+2k+l) 
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where i > 1 and only the first term on the right side of 3 is to be used when i = 1. It is 
convenient to set a^ = (-l)i + lbj, uj = (k+i)u-l, and i(i+2k+l) = 2ddj+i in equation 3, 
which leads to the formula 


d. . b. . = -u.b. + 
l+l l+l i i 


i-i 

E 

1 
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(4) 


From this result, it will be shown that b^ is a polynomial in u of degree i. We may 
take the first coefficient b* as an arbitrary but fixed nonzero constant. Then d 2 b 2 = 1 - 
u(k+l)b 1} a polynomial in u of degree one with the root u = Similarly, the first 

term on the right side of equation 4, in each case, is one degree higher in u than b^. 
For each root u of this polynomial, equation 4 holds for i = 1, 2, . . . , i-i and the 


equation 



(5) 


holds. For the system of equations 4 and 5 to hold simultaneously and yield non-trivial 
solutions, it is necessary and sufficient that the determinant of the coefficients vanish, 
namely 
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where we have set ej = — jh- . The determinant Pj in equation 6 is a polynomial in u of 
degree i and each root of b| + j = O is a root of P^= O, and conversely. Hence, the two 
polynomials differ only by a constant multiple. Comparing the leading coefficients, we 
readily deduce the relation 


P. = d 2 d 3 ... d. . b. .. 
l 1 6 l+l i+i 


(7) 


By a careful analysis of the determinant P i? it is found that the coefficients of u 1 , u i_1 > 
u 1-2 , and u* -2 are respectively 
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It follows that the sum of the roots of Pj is 2'7k+J) ’ ^he sum ^he P r °ducts of the 


roots taken two at a time is 
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which simplifies to 


i(i-l) 

2d(k+i) 


- Z 
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(k+fn) (k+j) 


the sum of the squares of the roots is 

i (i-1) 
d (k+i) 
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Z 


(k+j)' 


since ^ xj - ( ^ xj)" - 2 Y x m Xj. These considerations provide precise information 

pm i 

concerning the roots. The roots deviate about the values and the deviations have 

k+j 

a net sum of zero so that some are positive and some are negative (the positive devi- 
ations might be accounted for mostly because of pairs of conjugate imaginary roots) . 

The average deviation of the squares of the roots from the values 


is - 


i-i 


(k+j) 2 d(k+i) ’ 

which definitely suggests the presence of imaginary roots in many cases for i > 1. When 
i is increased by 1, a new root is introduced and the net effect on the sum of the roots is 

as though the roots of P. +1 are k+ . 1 +1 and the roots of Pj. However, each real root of 
p i is greater than — as is verified by an analysis of the determinant Pj, which has a 


positive value when u < 


k+i 


This may be shown by substituting such a value in Pj 


for u and successively eliminating the -dj, leaving only positive values down the main 
diagonal and zeroes everywhere below it; or, in case u =— , one can eliminate all 

K i 1 

the -dj except -d^ and, as before, Pj has a positive value. Since the eigenvalues are 
given by u 2 and £ u 2 = - + £ 1 


, we have this condition on the po- 


(k+j) 2 

tential eigenvalues. The sum of the products of the roots taken three at a time is 


Z 1 + L_ y 2 j ( j+1) (j+2k+l) ( j+2k+2) 

p> j> m (k+m) (k+j) (k+p) 8d 2 ^ ( j+k) (j+k+i) (j+k+2) ’ 

the second summation having the value i 1 ~ 2 ^ 1+1 > — when k=Q) so that the net deviation 

in this case is — . Considering all of these facts, there is an indication that 

for d large, the deviations are small and vice versa, that there are a few large positive 
deviations in the most recently introduced roots and in the case of conjugate imaginary 
roots, and that the less recently introduced roots (which are the largest) tend to stabi- 
lize slightly below a value ■ , V at the square root of a potential eigenvalue. One might 
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take u 2 = 


1 


(k+j ) 2 

of the squares of the roots from the values 
1 

~ d 


- as an approximation to an eigenvalue since the average deviation 

* is 7 — , which approaches 


(k+j ) 2 ~~ d(k+i) 

as i becomes infinite. This suggests the inequality d > (k+j ) 2 as a first approximate 


limit to the number of eigenvalues where d is finite. This will be reconsidered later. 


The coefficient 2dxe of one term in y or equation 2 is approaching zero at a 
weakly varying rate for extremely large values of x, when compared to the other coef- 
ficients. It might be expected that a reasonable approximation would be obtained by 
treating e " x as a constant, leading to the equation 


-c . 


i +1 


a i+l = [ (k+i) u " e 4 


( 8 ) 


and the eigenvalues 

e~ Cn 
U k+n 

where the notation c n indicates that the constant varies with n, which takes into account 
the fact that 8 is only an approximation to equation 3. Through the first two terms, when 

e n is expended into a Taylor series, this agrees with the discussion for k=0 of Ecker 

and Weizel provided we take c n equal to ^ (k+n) 2 . Thus we have as a second estimate 
to the eigenvalues 

i _ < k+n > 2 

u = — e d (9) 

k+n 


To check these approximations further, the determinant equation 6 was solved 
numerically for i = 1, 2, . . . , 10. In addition, the values so obtained were further re- 
fined by evaluations of the determinants P 2 q, P 40 J an< ^ ^ 41 * Finally a routine was devised 
based on the recursion formula 4 and the calculations made including the cases i = 1, 

2, . . . , 75. The latter routine was especially efficient, giving information in 75 cases 
instead of one case and in about one-fifth of the time taken for the one case. 


IE. CONCLUDING REMARKS 

On the basis of the above and other considerations, the following is presented as 
the best estimate of the values of the eigenvalues 

1 k+n ^ (k+n) 3 - 2k(n+l) 

U ~ k+n " ~ d + 4d 2 ' (10) 
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The values u of equation 10 are thought to be a little more than the true values, leading 
to the implication that the number of eigenvalues are limited by the following inequality 


2d s (k+n) 2 + | 2k(n+l) (k+n) I z , 


'll) 


where n is to take positive integral values and satisfy the inequality. The number of 
eigenvalues is not greater than the maximum such n (this assumes the correctness of 
the previous remark) . 

The numerical results obtained as explained above are summarized in the following 
tabular data. The values taken for d have been used on the basis of numerical conven- 
ience in machine computation. The values themselves are not as important as the gener- 
al magnitude of this parameter, which may take on any non-negative real value. The 
parameter k, on the other hand, must be a non-negative integer. The tabular data could 
be expanded on the basis of the techniques already developed. In fact, the predictor 
equation 10 may be sufficiently accurate for most cases, the truncation error apparently 

_ m o 

being of the order of d (this remark is based on theoretical reasons as well as on 
the tabular data below) . 

d = 1000 d = 100 


k = 0 

k = 1 

k = 2 

O 

II 

11 

k = 2 

. 99900025 



. 99002475 



.49800199 

. 49800099 


.48019431 

.48009421 


. 33034004 

. 3303384 

. 33033554 

.30396993 

.30381937 

. 30351756 

.24601583 

.24601383 

.24600983 

.21145270 

.21125064 

.21084493 

. 19503074 

. 19502823 

. 19502323 

. 15271597 

. 15246021 

. 15194527 

. 16071941 

. 16071640 

. 16071040 

. 12454326 

. 12285380 

. 12136559 

. 13594022 

. 13593672 

. 13592972 




. 11712289 

.11711888 

. 11711085 




. 10228432 

. 10227976 

. 10227081 


d = 10 


. 090234950 

. 090229537 

.090219821 




.080219801 

.080214067 

. 080201641 

O 

II 

k = 1 

<N 

II 

.071725530 

. 071697400 

. 071703049 

.90228380 



.064500631 

. 064424424 

. 064505604 

.40339323 

. 39722151 

. 37764286 


In general, as i increases in Pj_, the roots for k and n large make their initial ap- 
pearance. The tabular data indicates a definite trend in agreement with inequality 11 
and equation 10 agrees with the data in a remarkable way, as one may verify by evalu- 
ation of u in equation 10. The estimate due to Ecker and Weizel is equivalent to 

3 3 

i - 5 + -‘V- - . . . , if we take x.. to be "r - , which gives the most favorable choice in 
n d 6d 4 u cl 
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general. It may be seen from the tabular data that by our choice of Xq, the results are 
good but do not give an estimate as good as equation 10. 

For each such eigenvalue, the eigenfunction F n is given by the equation 

F„ - e -udx l (-l) J+1 b x k+ S, (12) 

3=1 

where bj is an arbitrary constant not zero, bj + j = Pj (d 2 d 3 . . . dj + ^) ^ , d^ + ^ = 

j(j+2k+l) , is the determinant in equation 6, u is determined by equation 10 (approxi- 
mately) , and n satisfies the restriction 11. 
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